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© Continuous core samples are obtained from 
along the length of a wellbore. Discrete samples are 
obtained from the continuous core samples at se- 
lected intervals so as to be representative of the 
formation lithologies transversed by the wellbore 
along its iength. The samples can then be processed 
at the weilsite with a transportable laboratory so as 
to obtain measures of a plurality of selected physical 
properties of the subterranean formations represen- 
tative of their in-situ condition. An inverse log of 
each physical property as a function of depth can be 
prepared with which explorationists can evaluate 
geophysical survey data for oil and gas potential. 
Additionally, inverse logs of one or more of the 
physical properties can be used to perform basis 
analyses. When used for basis analysis, the physical 
properties of each inverse log can be used to cali- 
brate geophysical survey data of the basin at one or 
more wellbore sites within the basin from which the 
inverse logs were obtained. 
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A METHOD AND SYSTEM F 

The present invention relates generally to a 
method and system for geophysical exploration 
and more particularly to a method and system for 
obtaining selected physical properties from sam- 
ples of the earth's subterranean formations sur- 
rounding a wellbore and for obtaining inverse logs 
of the selected physical properties along the length 
of the wellbore for evaluating the oil and gas bear- 
ing potential of the earth's subterranean formations. 

in the continuing search for oil and gas. ex- 
plorationists have developed a wide array of geo- 
physical exploration methods for imaging the 
earth's subterranean formation structures and dis- 
playing those images as well as providing es- 
timates of selected formation characteristics such 
as lithology. porosity, density, and velocity of prop- 
agation of seismic energy in the earth's formations. 
Exemplary of such geophysical exploration meth- 
ods seismic surveys, gravity surveys, and magnetic 
surveys. The principal geophysical method for oil 
and gas exploration within the earth's subterranean 
formation is - seismic exploration. Simplistically, 
seismic exploration comprises imparting seismic 
energy into the earth and recording the earth's 
response thereto. The earth response can include 
both refracted and reflected seismic energy. As the 
seismic energy propagates through the earth, it is 
controlled by the physical properties of the rocks 
and the fluids contained therein. Unfortunately, ex- 
plorationists cannot directly measure these phys- 
ical properties which constrain the propagation of 
seismic energy from the recorded seismic survey 
data. Rather, seismic exploration techniques mere- 
ly provide a measure of the time it takes the 
seismic energy to traverse some path. In fact, it is 
often times impossible to precisely ascertain the 
path taken by the seismic energyi In view of such 
uncertainties, how is it that explorationists are able 
to estimate formation lithology, fluid content or 
even simple, reconstruct formation geometries 
through which the seismic energy propagates with, 
only a measure of time? Quite simply, the answer 
is explorationists cannot without the aid of addi- 
tional information such as estimates of the veloci- 
ties of propagation of seismic energy In rock and 
rock density as well as other rock and fluid phys- 
ical properties. Methods for estimating the velocity 
of propagation and a density as a function of depth 
are generally known to those skilled, in the art; 
however, velocity estimations are the weakest ele- 
ment in the inversion of seismic time data. 

Often times, wells which have been drilled 
based on geophysical survey data are logged us- 
ing one or more of the well logging techniques 
such as acoustic logging, resistivity logging, neu- 
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tron logging, etc. Results of such well logging 
techniques can provide the exploration! st with se- 
lected physical properties generally over regions of 
interest along the wellbore length. With such well 

5 logging techniques, explorationists are able to ob- 
tain better estimates of the formation physical 
properties for evaluating geophysical survey data. 
However, it must be recognized that such well 
logging techniques still only provide an indirect 

w measure of the formation physical properties of 
interest and are generally limited to only selected 
regions of interest within a wellbore. 

Occasionally, core samples are obtained from 
limited regions of interest within a wellbore for 

75 transport to a central laboratory so that physical 
properties of the formations adjacent the borehole, 
such as fluid content of the sample, (U.S. Patent 
2,458,093 M. Muskat, et al., January 4. 1949), oil 
content of the sample (U.S. Patent 2,531.083, A. L. 

20 Smith, November 21, 1950), and electrical resistiv- 
ity of the sample (U.S. Patent 2.613.250, H. L. 
Bilhartz, et al. a October 7, 1952), can be obtained. 
Unfortunately, such core samples generally repre- 
sent only a small fraction of the formational 

25 lithologies traversed by the wellbore along its 
. length. Moreover, the delay associated with trans- 
porting the samples to the central laboratory can 
preclude the measurement of properties of certain 
friable lithologies or result in seriously degraded 

30 measurements of certain properties. During the 
drilling process, chips which have broken away 
from the subsurface formations can be captured 
from the drilling fluid at the surface for analysis. 
Typically, explorationists have identified the 

35 lithologies of the subsurface formations from such 
chips. However, determining the exact depth in the 
wellbore from which the chips originated Is gen- 
erally not possible due to the very nature of their 
transport to the surface. Consequently, core sam- 

40 pies and chip samples fail to provide either a 
regular sampling of borehole formation lithologies 
along the length of the wellbore or fail to accurately 
determine the depth of origin of such samples so 
that logs of formation characteristics, as a function 

4$ of depth over generally the entire wellbore length, 
cannot be produced* Moreover, the very nature of 
the chip samples can often times preclude measur- 
ing certain formation characteristics (for example, 
measurement of rock velocity or compressive 

so strength). 

In spite of current geophysical exploration 
methods, explorationists are not provided with di- 
rect measures of the physical properties of the 
formations surrounding a wellbore along its length. 
Rather, attempts to evaluate and~thterpret geo- 
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physical survey data must generally depend on 
estimates of formation physical properties or on 
indirect measures of the formation physical prop- 
erties from only limited regions of interest in the 
earth's subsurface. Consequently, explorationists 
are unable to accurately calibrate the vertical extent 
of their geophysical survey data with direct mea- 
sures of the earth's physical properties so as to 
better evaluate the oil and gas bearing potential of 
the earth's subterranean formations. 

In order to overcome the difficulties noted 
above, the present invention provides a novel 
method and system of obtaining selected physical 
properties from samples of the earth's subterra- 
nean formations surrounding a weilbore and devel- 
oping inverse logs of the selected physical prop- 
erties along the length, of the weilbore. Expioration- 
ists can employ such inverse logs of the selected 
physical properties to perform basin analysis 
whereby the physical properties obtained at one or 
more wellbores can be used to calibrate geophysi- 
cal survey data at each weilbore location and then 
extrapolate such physical properties to other loca- 
tions in the basin with the geophysical survey data 
so as to evaluate the oil and gas bearing potential 
of the earth's subterranean formations. 

The present invention relates generally to a. 
method and system for geophysical exploration 
and more particularly to a method and system for 
obtaining selected physical properties from sam- 
ples of the earth's subterranean formations sur- 
rounding a weilbore and for obtaining inverse logs 
of the selected physical properties along the length 
of the weilbore so as to aid the geophysicist in 
evaluating the oil and gas bearing, potential of the 
earth's subterranean formations. 

In one embodiment of the invention, discrete 
samples can be obtained at selected intervals 
along the length of a weilbore so as to be repre- 
sentative of the formation lithplogies traversed by 
the weilbore along its length. The discrete samples 
can then be processed at the wellsite with a trans- 
portable laboratory so as to obtain measures of a 
plurality of selected physical properties of the sub- 
terranean formations representative of their in-situ 
condition. An inverse log of each physical property 
as a function of depth can be prepared with which 
explorationists can evaluate geophysical survey 
data for oil and gas potential. In more detail, mea- 
sures of the subterranean formation's compres- 
sional wave velocity, shear wave velocity, grain 
density, saturated density, dry density, magnetic 
susceptibility, compressive strength, and static 
elastic moduli can be determined from which one 
can also calculate the formation's, acoustic imped- 
ance and porosity. 

Inverse logs of one or more of the physical 
properties can be used to perform basin analyses. 



When used for basin analysis, the physical prop- 
erties of each inverse log can be used to calibrate 
geophysical survey data of the basin at one or 
more weilbore sites within the basin from which the 

5 inverse logs were obtained. Having thus calibrated 
the* geophysical survey data, explorationists can 
then extrapolate the physical properties to other 
locations in the basin remote to the weilbore site 
with the geophysical survey data to evaluate the oil 

w and gas potential of the earth's subterranean for- 
mations at such remote locations. 

In another embodiment of the invention, a 
transportable laboratory is provided which is espe- 
cially adapted for processing, at a wellsite. large 

fs numbers of discrete samples obtained along the 
length of a weilbore in a minimum amount of time 
so as to obtain measures of the formations physical 
properties which more properly reflect their in-situ 
condition. Timely measurement of the physical 

20- properties of certain friable lithologies (such as 
shales) is essential if physical properties are to be 
obtained at all. Additionally, the in-situ measures of 
certain physical properties can be seriousJy deg- 
raded by a deJay in measurement e.g.. seismic 

25 energy velocities of propagation in the rock, and 
porosity. Surprisingly, such properties have been 
found to vary over a wide range of values simply 
due to a delay in the time of their measurement. 
The transportable laboratory comprises a plu- 

30 rality of highly automated measurement stations 
each connected to a central processing unit (CPU). 
By labeling each sample with a unique machine 
readable sample code and each measurement sta- 
tion with a unique machine readable station' code, 

35 reader means can be provided for reading the 
sample code of each sample and the station code 
of a selected measurement station into the CPU. 
Upon receipt of valid sample and station codes, the. 
CPU Initiates the measurement of the sample at 

40 the selected measurement station and stores the 
measurements \n a data file associated with the 
unique sample code. The CPU can also calculate 
additional physical properties from selected mea- 
sured physical properties. Thus, the transportable 

45 laboratory becomes In effect a "paperless" labora- 
tory since manual recordation of a sample's Iden- 
tity and it's measured and calculated physical 
properties has been eliminated. The paperless lab- 
oratory greatiy reduces errors generally associated 

so with manually recording such Information. In addi- 
tion to Initiating such measurements and recor- 
dation, the CPU can initiate a prescribed sequence 
of measurements at a given measurement station 
as well as collect, collate and evaluate the physical 

55 properties measured to provide inverse well logs of 
the measured properties as a function of depth. 
Additionally, the CPU can calculate additional phys- 
ical properties from the physical properties mea- 
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sured as well as reinterpret geophysical survey 
data and models at a remote site. 

The. present invention will be better understood 
with reference to the following drawings and de- 
tailed description- In the drawings:- 

Figure 1 is a schematic representation of a 
transportable laboratory according to the present 
invention; 

Figure 2 is a generalized flow diagram for 
the operation of each measurement station in the 
transportable laboratory; 

Figures 3 (a. b, c) depict inverse logs of 
velocity measured from samples obtained at se- 
lected intervals along three separate wellbdres; 

Figures 4 (a, b, c) depict inverse logs of 
porosity measured from samples obtained at se- 
lected intervals along three separate wellbores; 

Rgures 5 (a, b, c, d, e) depict inverse logs of 
formation properties density, compressional veloc- 
ity (Vp), shear velocity (Vj) Vp/V s ratio. Impedance 
and magnetic susceptibility from samples obtained 
at generally ten foot intervals along the length of 
the welibore; and 

Rgures 6 (a, b, c, d, e) depict smoothed 
inverse logs of density, compressional velocity (v p ), 
shear velocity (v s ), the ratio Vp/v 8l and formation 
impedance, respectively, while Figure 6 (f) repre- 
sents synthetic seismic traces obtained by convolv- 
ing a source wavelet with the inverse log of forma- 
tion impedance shown in Rgure 5 (e); • 

Rgures 7 (a, b) Show a section of a seismic 
survey in Rg. 7a, and Rg. 7b depicts an inverse 
log of formation impedance obtained from samples 
collected along the length of the welibore indicated 
by the derrick structure in Rg. 7a. 

The present invention provides a novel method 
and system for geophysical exploration. 

The following discussion is provided in order to 
more completely understand the present invention. 
In the continuing search for oil and gas, ex- 
plorationists employ a variety of geophysical and 
geological survey techniques to identify regions 
having potential oil and gas formations. By their 
very nature, such geophysical and geological sur- 
vey techniques generally provide only indirect es- 
timates of selected subsurface formation properties 
which the explorationist uses to interpret potential 
oil and. gas formations. Moreover, such survey 
techniques generally depend on at least estimates 
of selected physical properties of the earth's sub- 
terranean formations before any meaningful inter- 
pretation can be made of either the formations 1 
structure or their oil and gas bearing potential 

By way of example, explorationists often times 
employ the time it takes seismic energy to travel 
from the earth's surface to a reflecting horizon and 
back again to infer structure, lithology, porosity 
fluid content as well as the path traversed by the 



seismic energy. Similarly, magnetic surveys can 
provide estimates of formation structure. However, 
geophysical and geological survey techniques do 
not generally directly measure the properties of 

5 interest, but rather simply provide indirect mea- 
sures of other characteristics of the formation from 
which it is possible for the explorationist to es- 
timate the physical properties of interest. 

After a welibore is drilled, well logging tech- 

w niques (such as sonic logging, electric resistivity 
logging, gamma ray logging, gravity logging, neu- 
tron logging, etc.) can provide additional estimates 
of the subsurface formation properties of interest 
However, additional sources of error (such as 

/$ borehole rugosity, mudcaking, varying borehole di- 
mensions, variation of and borehole fluids) can be 
introduced by such logging techniques which may 
adversely affect the quality of the estimates of the 
formation physical properties. On occasion, well- 
20 bores have been cored to provide a few samples of 
selected formations over limited intervals adjacent 
the welibore from which direct measures can be 
obtained. Generally, such samples are obtained at 
a wellsite and then transported to a central facility 

25 for measurement some time thereafter. However, 
Applicants have found that certain formation prop- 
erties {such as velocity and porosity) can be sub- 
stantially altered if they are not measured shortly 
after removal from the welibore, while the physical 

30 properties for certain friable lithologies are simply 
not obtainable at all if not done within a few hours 
after their removal from the welibore. Moreover, 
even for those formation physical properties not 
affected by time, the physical properties measured 

35 are generally measured on only a handful or so of 
samples over a very short interval of the welibore. 
Consequently, the explorationist is simply not pro- 
vided with direct measures of formation physical 
properties for substantially the entire length of the 

40 welibore. Hence, the value of such information is 
limited to regions adjacent the borehole from which 
the samples came and the explorationist is left to 
evaluate and interpret the balance of his geophysi- 
cal survey data with only estimates of the subterra- 

45 nean formation physical properties. 

In one embodiment of the invention, core sam- 
ples can be obtained over substantially the entire 
length of a welibore; however, an important aspect 
of the present invention is that discrete samples of 

50. the earth's formations surrounding a welibore need 
only be taken at selected intervals along the length 
of the welibore. One criteria for establishing the 
sampling intervals is that discrete samples gen- 
erally be obtained for every formation lithology 

55 traversed by the welibore along its length. Since 
the physical properties of the formations surround- 
ing the welibore can also be used to evaluate and 
calibrate geophysical survey data_Je.g., seismic 
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survey data), the minimum resolution interval of the 
geophysical survey data being evaluated estab- 
lishes another criteria for the sampling interval. 
Typically seismic data has the finest resolution 
interval of any geophysical survey technique 5 
(approximately 100 ft). In practice, we have found 
that a sampling interval of one sample every 10 ft 
adequately addresses the criteria of lithology 
changes as well as the geophysical survey mini- 
mum resolution interval criteria so as to address w 
the problems of heterogeneous rather than homo- 
geneous earth formations so that the physical prop- 
erties of the samples can be correlated with geo- 
physical survey data without introducing sampling 
bias. ,5 

In another embodiment of the Invention, a 
transportable laboratory is provided so as to pro- 
cess samples obtained from the earth's formations 
surrounding a wellbore along its length. Since time- 
ly measurement of certain physical properties and 20 
of certain lithologies can be important the trans- 
portable laboratory is preferably configured for field 
operation in the vicinity of the wellbore. Generally, 
the transportable laboratory comprises a plurality of 
measurement stations, each adapted to measure at 2s 
least one physical property of each sample. Addi- 
tionally, the transportable laboratory includes a 
central processing unit for identifying samples, ini- 
tiating measurements on the samples at selected 
measurement stations, and storing the measured jo 
physical properties in a sample data file. 

Prior to discussing the operation of each mea- 
surement station of the transportable laboratory A, 
the basic sequence of steps for all operations will 
be discussed with reference to both Figures 1 and 35 
2. Every sample which enters the transportable 
laboratory A is given a unique machine readable 
sample code (which includes a we If identifier, date 
and depth). Preferably, such sample codes are 
included on a label adapted to be mounted on 40 
each sample. Every measurement station is as- 
signed a unique machine readable station code. 
Prefera bly, the sample codes and station codes 
are bar codes which can be read by a laser bar 
code reader into a central processing unit (CPU). 45 
As such, the transportable laboratory A is capable 
of processing a large number of samples through 
each of the various measurement stations employ- 
ing multiple, asynchronous flow paths, until each 
sample has been through each measurement sta- so 
tion. To initiate a measurement sequence at any 
measurement station/ the machine readable station 
code (preferably a bar code)^ assigned to a se- 
lected measurement station is read into the CPU 
using a laser bar code reader 70 at step 100. The 56 
sample to be measured at the selected measure- 
ment station then has its unique machine readable 
sample code (preferably a bar code) read Into the 



CPU with the laser bar code reader 70 at step T10. 
If the station code and the sample code are both 
valid codes, the CPU directs the selected measure- 
ment station at Step 120 to initiate a sequence of 
measurements unique to the selected measure- 
ment station and the set of measurements obtained 
are then stored in the CPU at step 130 in a sample 
data file associated with the sample's unique sam- 
ple code. 



I. SYSTEM 

Looking now to Figure 1, a schematic repre- 
sentation of the transportable laboratory A is 
shown. Although the transportabfe laboratory A is 
particularly suited to field needs, those skilled in 
the art will appreciate that the transportable labora- 
tory satisfies many of the needs of a central sta- 
tionary laboratory. Generally, the transportable lab- 
oratory A comprises a plurality of uniquely iden- 
tified measurement stations connected to a central 
processing unit (CPU) such as a Hewlett-Packard 
9000 series 500 computer and a laser bar code 
reader 70 as manufactured by Intermec: The CPU 
also includes a mass storage system (not shown) 
for collecting and storing the physical properties 
obtained from the various samples which are. pro- 
cessed by the transportable laboratory A. 

In particular, the transportable laboratory 5 A in- 
cludes an electronic dimension jig 10 for digitally 
measuring sample lengths and diameters. The di- 
mensioning jig 10 includes pairs of opposed linear 
voltage displacement transducers (LVOT) adapted 
to measure axial and lateral dimensions of each 
sample. The transportable laboratory A afso in- 
cludes an electronic balance 20 adapted to mea- 
sure the weight of each sample. Further, an elec- 
tronically driven pyqnometer 30 is provided to 
measure grain density of each sample; a triaxial 
loading apparatus 40 is provided to obtain mea- . 
sures of each sample's compressive strength, as 
well as provide an input for calculating Young's 
modulus and Poisson's ratio for each sample; an ' 
electronic magnetic susceptibility meter 50 is pro- 
vided to obtain a measure of the magnetic suscep- 
tibility for a given weight of sample; and a pair of 
electronically driven pressure velocity cells 60 are 
provided to measure the velocities of propagation 
of compression^ (P), horizontal shear (SH), and 
vertical shear (SV) waves in each sample as a 
function of pressure. A laser bar code reader 70 is 
provided to read both the unique sample codes 
assigned to each sample and the unique station 
codes assigned to each measurement station into 
the CPU such that the CPU can initiate a measur- 
ing sequence on a sample at a selected measure- 
ment station and store the results of such measure- 
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merit in a sample data file within the CPU memory. 
The CPU accumulates all the measurements from 
the various measurement stations for a particular 
sample in the sample data file. The individual sam- 
ple data files can be displayed at a printer or 5 
plotter 80 or a CRT screen 90. or be transmitted 
via satellite. to a central location with a satellite 
communication station 95. Additionally, each phys- 
ical property from the samples can be sequenced 
according to depth to produce inverse logs of se- w 
lected physical properties representing the forma- 
tions surrounding the wellbore along its length or 
used form cross plots of selected physical prop- 
erties for trend analysis. 

The transportable laboratory A of the present /5 
invention provides means for obtaining measure- 
ments of the following physical properties from 
discrete samples taken at selected intervals along 
the length of a wellbore: grain density, dry bulk 
density, saturated bulk density, magnetic suscepti- 20 
bility, compressive strength, elastic moduli, and 
seismic velocities (P. SH. SV) versus pressure. 
With these direct measurements one can calculate 
many other physical properties of interest For ex- 
ample/ the difference between the grain density 25 
and the dry bulk density yields a porosity estimate; 
likewise, the difference between the saturated and 
dry bulk densities yields another estimate of poros- 
ity. Further, combinations of saturated density and 
seismic velocities ' yield dynamic elastic moduli, '30 
acoustic impedance and reflection coefficients. 
Piecewise integration of the seismic velocity over 
the depth interval sampled can yield an estimate of 
the two-way traveltime, -which can be used to cali- 
brate the time observed on a seismic section to a 35 
true depth. Additionally, synthetic seismic traces 
can be produced by convolving an inverse log of 
formation impedance with a wavelet having known 
characteristics, for the calibration of seismic survey 
data. A more complete discussion below will detail to 
uses of the measured physical properties and the 
inverse logs derived therefrom. 



II. OPERATION 45 

Each discrete sample obtained at selected in- 
tervals along the length of a wellbore is assigned a 
unique machine readable sample code, preferably 
a bar code, which can be mounted with the sample 50 
or attached with its sample holder. All measure- 
ment sequences begin by entering the measure- 
ment station code and sample code via the bar 
code reader 70 into the CPU. Hereafter, the step of 
reading the station code and sample code will be 55 
understood in the operation of each measurement 
station. In order to carry out the desired measure- 
ments at the various measurement stations shown 



in Figure 1. three separate samples~_are cut from 
each discrete sample. The three samples obtained 
from each discrete sample are each generally right 
circular cylinders having lengths of approximately 
13mm. 38m and 70mm (0.5", 1.5" and 2.75"). With 
reference to Figures 1 and 2, generally the detailed 
steps to obtain at least one physical property of 
each sample at each measurement station will be 
discussed. In order to. simplify the following discus- 
sion, the operation of the transportable laboratory A 
will be made with respect to the operations which 
occur on each separately sized sample 

The 13mm sample will be referred to hereafter 
as the disc sample. The initial measurements made 
on each disc sample are its length and diameter at 
10 and its saturated weight at 20. The disc sample 
is dried and then placed In a vacuum oven for an 
additional period. The dried disc samples are then 
weighed again at 20 to obtain their dry weight. The 
dry weight can be used by the CPU in calculating 
the dry bulk density of each sample. About 20 gms 
of the disc sample are subsequently crushed to a 
particulate size of 6mm (0.25") or less. A portion of 
the fragmented disc sample is then placed "in a 
pycnometer cup and reheated in a vacuum oven. 
The pycnometer 30 is used to obtain a measure of 
the disc sample's grain density. The pycnometer 
30 measures a volume of space occupied by the 
fragmented sample and pycnometercup using he- 
lium gas. Because the measurements by the pyc- 
nometer 30 takes the longest time of any measure- 
ment station, the system includes multiple pyc- 
nometers 30 which are each interfaced to an IBM 
AT computer (IBM) (not shown) so as to permit 
automated operation of each pycnometer which in 
turn is connected with the CPU. 

After receiving a valid sample code and pyc- 
nometer station code, the CPU signals the selected 
pycnometer to begin its measurements. The frag- 
mented disc sample is inserted into the pycno- 
meter's sample chamber. Upon completion of the 
measurement cycle, the IBM reads a measured 
volume and divides by the recorded sample weight 
to arrive at a grain density for the fragmented disc 
sample. This value along with the sample code and 
pycnometer station code are sent back to the CPU 
for storage and further processing. 

The fragmented disc sample can also be used 
for magnetic susceptibility measurement at 50. 
Here the fragmented disc sample is poured into a 
polyethylene vial of known weight and reweighed at 
20 and its specific magnetic susceptibility is mea- 
sured using a Bartington magnetic susceptibility 
meter at 50. When the magnetic susceptibility 
measurement has been completed, an average of 
five readings can be sent to the CPU memory. 

The 38mm sample can be used to measure 
the seismic velocities of the discrete samples at 
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pressures appropriate for their in situ depth. The 
magnitude of seismic velocity in rock is sensitive 
the concentrations of cracks and pores. Since 
cracks exhibit a strong, pressure dependence, seis- 
mic rock velocities can exhibit a strong pressure 
dependence. Therefore, to make meaningful ve- 
locity measurements on rock, one must subject it 
to the pressure and saturation it experiences in 
situ. This is done most effectively by making ve- 
locity measurements on the 38mm samples at a 
plurality of pressures which bracket the in-situ 
pressure appropriate for each discrete sample. 

The flat faces of the 38mm samples are first 
polished to remove cutting marks and to make the 
faces more parallel for dimension measurements at 
10 and velocity measurements at 60. A bar code 
label can be affixed to the perimeter of the 38mm 
sample with mylar tape. One side of the 38mm 
sample is marked to indicate a preferred bedding 
direction for alignment with respect to the shear 
wave polarizations. The 38mm samples are kept in 
a water bath whiie awaiting velocity measurement 
so as to maintain their fluid saturation. The 38mm 
samples are typically processed in batches since 
each pressure velocity,- cell 60 has positions for 
multiple samples. 

The 38mm sample processing begins by mea- 
suring its dimensions at measuring station 10. The 
38mm sample weight can then be obtained at 
measurement station 20. Having the 38mm sample 
dimensions and weight, one is ready to measure 
the velocities of the compressional (P) and shear 
(S) waves in the 38mm samples as a function of 
pressure using a pulse transmission technique. 

Each 38mm sample is inserted between a set 
of end caps. The relative orientation of each 38mm 
sample with the end cap provides a mechanism for 
determining the relative anisotropic orientation of 
the sample. Piezoe lectric transducers assemblies 
are mounted within , both end caps. The 38mm 
sample and end cap assembly are then Installed in 
the pressure velocity cell 60 and electrical connec- 
tors are attached. These connectors, permit the 
transmitted and detected signals to be applied and 
received from the appropriate transducers. The op- 
erator selects a pressure velocity cell by reading 
the appropriate station code and sample code via 
the bar code readers. Once the CPU decides there 
is a valid sample code, a control program begins 
filling and pressurizing the pressure velocity cell 
60. An acquisition program is activated by the GPU 
when equilibrium at programmed pressure set 
points Is achieved. Five time series are collected 
for each sample at each pressure point. When the 
measurements at the pressure velocity cell 60 are 
complete, the samples are removed and rewelghed 
at 20. The new weight is the saturated weight. 
Saturated weight plus sample volume obtained 



from sample dimensions yields a saturated density. 

With the measurement of. the compressional 
(P) and shear (S) wave velocities of each 38mm 
sample, the CPU can calculate a measure of the 

s 38mm. samples anisotropy, the magnitude of the" 
difference between the fast and slow shear veloci- 
ties, and saturated density. When these measure- 
ments are coupled with disc sample measure- 
ments, one can estimate sample porosity i.e., dif- 

10 ference between saturated and dry density provide 
a saturated porosity estimate. 

The remaining 70mm sample (hereafter plug) 
is plugged parallel to its long axis with a. 25mm (1 ") 
diameter hollow core bit. The ends are trimmed 

/s parallel and flat with a slow speed precision cutoff 
saw. The plug sample is used in the triaxial loading 
apparatus 40 to determine its compressive 
strength, Young's modulus, and Poisson's ratio. 
The plug sample length and diameter are mea- 

20 sured at 10. A radial LVDT is positioned at the 
midpoint of the plug sample and an axial LVDT is 
mounted on one end of the plug sample and exten- 
sion rod on the other. The plug sample can then be 
inserted Into the triaxial loading apparatus 40. The 

25 CPU first sets a confining pressure within the tri- 
axial loading apparatus 40 to a value appropriate 
for the sample's depth within the wellbore and then 
raises the pore pressure to a fixed value to assure, 
sample saturation. The CPU then initiates axial 
• 30 loading of the plug sample at a programmed rate, 
nominally 10-6 sec-1, until failure of the plug, sam- 
ple. Axial and radial displacements from the 
LVTS's are constantly recorded by the CPU during 
loading. These values can be normalized by the 

35 original dimensions to provide axial and radial 
strains. The CPU can also determine the slope of 
the axial strain versus axial stress to yield a value 
for Young's modulus. The incremental ratio of ra- 
dial to axial strain can provide values of Poisson-s 

40 ratio. The maximum stress (force/sample area): 
reached during the test is the compressive strength 
of the plug sample. 

Separate operating programs allow the CPU to 
process a plurality of samples through multiple, 

45 asynchronous flow paths until each sample has 
gone through each of the various measurement 
stations. These programs run in there own virtual 
machines within the CPU. 

The most complicated program is MAIN. This 

so program directly interfaces to the bar code readers 
70, magnetic susceptibility meter 50, balance 20, 
pycnometer, digital dimension jig 10 and memory 
storage devices (not shown) and indirectly with the 
pressure velocity cell 60 and the triaxial loading 

55 apparatus 40. Ail bookkeeping and sample data file 
entries are done through the program MAIN. When 
a sample code is recognized by the CPU, It. is 
assigned an external sequence number which con- 
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sists Of year, Julian day and daily run number. An 
active discrete sample code list is maintained by 
MAIN. The sample code list includes filenames and 
mass storage locations, The filenames correspond 
to the measurement results from selected mea- 
surement stations. When a sample code is entered 
by the laser code reader 70. MAIN first checks the 
active sample code list: if the sample code exists, 
it will proceed to the proper measurement subrou- 
tines. If the sample code is not found on the active 
sample code list, a set of data files are created and 
the sample code is added to the active sample 
code list and measurements can then be carried 
out. Whenever the active sample code list is up- 
dated or a new measurement is performed, the 
sample codes and status can be displayed on the 
CRT 90. 

There are a few measurements which are more 
involved than weighing or measuring sample di- 
mensions such as velocity and compressive mea- 
surements. For such measurements, additional pro- 
grams which vary extrinsic variables and programs 
to measure intrinsic properties of the samples as a 
function of the variation in extrinsic variables are 
provided. 

III. FORMATION EVALUATION 

With the direct measures of a plurality of se- 
lected physical properties of the formation sur- 
rounding one or more wellbores along their length, 
exploratfonists are now in a position for the first 
time to more intelligently use geophysical survey 
data to evaluate the oil and gas bearing potential of 
the earth's subterranean formation both during the 
well drilling process at the welisite as well as after 
completion of the well at a remote facility. Looking 
first to Figures 3, 4, 5, and 6, inverse logs of 
measured and calculated physical properties ob- 
tained from samples of the earth's formations at 
selected intervals along the length of a wellbore are 
shown. In particular, Figures 3 (a, b and c) each 
represent sent inverse logs of the compressional 
wave velocity (V p ) in samples obtained at selected 
intervals along the length of three separate wells. 
The velocities depicted In Figure 3 represent the 
velocity in a particular sample unlike previous ve- 
locity estimates which are really no more than an 
averaged or integrated velocity. The importance of 
such distinction is that acoustic impedance 
changes in the earth which provide an indication of 
the nature of the reflection and refraction of seismic 
energy depend not on averaged or integrated ve- 
locities; but rather, on the interval velocity of seis- 
mic propagation in each formation lithology. 

Figures 4 (a, b, and c) depict Inverse logs of 
formation porosity measured from samples ob- 



tained at selected intervals- along the length of 
three separate wellbores. Additionally, Figures 5 (a, 
b, c, d, e and f) depict inverse logs of formation 
density, compressional wave velocity (V p ), shear 

5 wave velocity (V s ), the ratio Vp,V„ formation imped- 
ance and magnetic susceptibility, respectively. The 
sample interval for the inverse logs of Figure 5 is 
much finer than the minimum resolution thickness 
of conventional geophysical survey data. Conse- 

io quently, such oversampling tends to produce in- 
verse logs which have a very "spiky" appearance 
as depicted in Figure 5. However, when the inverse 
logs of Figure 5 are smoothed with a smoothing 
operator having a wave length generally equal to 

/5 the minimum resolution thickness of the geophysi- 
cal survey data with which it is to be compared, the 
inverse logs take on. an appearance similar to those 
in Figures 6 (a-e). In fact, the inverse logs of Figure 
6 (a-e) were smoothed with a smoothing operator 

20 having a minimum resolution thickness appropriate . 
for seismic survey data. The synthetic traces of 
Figure 6 (f) were produced by convolving a wavelet 
signal, having characteristics of a seismic source 
with the smoothed inverse impe dance log of Fig- 

25 ure 5 (e) and correspond well with actual seismic 
survey data. 

The importance of having such inverse logs of 
formation properties along the length of a well both 
during the drilling process and upon completion of 

so the drilling process can be very important By way 
of example with the inverse logs of formation ve- 
locity to the point of current drilling within the 
wellbore. the original geophysical survey data can 
be reprocessed to forecast depth to target Addi- 

35 tionally, Inverse logs of formation compressive 
strength obtained during the drilling process can 
provide the drillers with information important for 
controlling drilling parameters. 

With inverse logs of formation physical prop- 

40 erties obtained at selected intervals along the well- 
bore, the geophysical survey data can be cali- 
brated to the physical properties represented in the 
inverse log. In particular,. Figures 7a and 7b repre- 
sent a seismic time section and an inverse log of 

4s formation Impedance measured from samples ob- 
tained at selected intervals along the length of a 
wellbore. In particular, the inverse impedance log . 
of Figure 7b was obtained at a well site indicated . 
by the derrick structure in Figure 7a, With such 

so comparison of the inverse log of impedance and 
the seismic section one can calibrate the physical 
properties of the formation at the location of the 
wellbore site. Having calibrated the geophysical 
survey data (in this case the seismic survey data) 

55 at the wellbore site, geophysicists can then use the 
calibrated geophysical survey data to extrapolate 
the physical properties to remote locations in^ a 
basin, 



8 



15 



EP 0 367 495 A2 



16 



Thus, the present invention provides a novel 
method and system of obtaining selected physical 
properties from samples of the earth's subterra- 
nean formations surrounding a wellbore and devel- 
oping inverse logs of the selected physical prop- 
erties along the length of the wellbore. Exploration- 
ists can employ such inverse logs of the selected 
physical properties to perform basin analysis 
whereby the physical properties obtained at one or 
more wellbores can be used to calibrate geophysi- 
cal survey data at each wellbore location and then 
extrapolate such physical properties to other loca- 
tions in the basin with the geophysical survey data 
so as to evaluate the oil- and gas-bearing potential 
of the earth's subterranean formations. 

It should thus be understood that the invention 
is not limited to any specific embodiment set forth 
herein as variations are readily apparent and thus 
the invention is to be given the broadest possible 
interpretation within the terms of the following 
claims. 



Claims 

1 A geophysical system for processing a plu- 
rality of samples of the earth's subterranean forma- 
tions so as to obtain a plurality of physical prop- 
erties for each sample representative of their in-situ 
condition, comprising: 

(a) a unique machine readable sample code 
label with each sample; 

(b) a plurality of measurement stations for 
measuring a plurality of selected physical prop- 
erties, wherein each measurement station is adapt- 
ed to measure at least one physical property for 
each sample and each measurement station in- 
cludes a unique machine readable station code; 

(c) reader means for reading the sample 
codes and station codes for a sample at a selected 
measurement station; and 

<d) central processing means for initiating 
the measurement of at least one physical property 
of a sample at a selected measurement station 
upon receiving the station code and the sample 
code from the the reader means as well as for 
storing the measurements of the sample. 

2. The geophysical system of Claim 1 further 
including: 

display means for displaying inverse togs of se^ 
lected physical properties as a function of depth in 
a wellbore. 

3. The geophysical system of Claim 1 further 
including: 

means for transmitting the measures of physical 
properties to a central location. 

4. The geophysical system of Claim 1 wherein 
the measurement stations include means for mea- 



suring formation density, compressional wave ve- 
locity, shear wave velocity, compressive strength, 
porosity and magnetic susceptibility. 

5. A method for processing a plurality of sam- 
5 pies of the earth's subterranean formations at a 

plurality of measurement stations so as to- obtain a 
plurality of physical properties for each sample 
representative of their in-situ condition, comprising 
the steps of: 

/o (a) labeling each sample with a unique ma- 

chine readable sample code; 

(b) assigning each measurement station a 
unique machine readable station code; 

(c) machine reading the sample code and 
is station code of a sample at a selected measure- 
ment station into a central processing unit; and 

(d) having the central processing unit re- 
sponsive to receipt of the sample code and the 
station code initiating a measurement of at least 

so one physical property of the sample at the selected 
measurement station and storing the measure of 
the physical property. 

6. The method of Claim 5 further include the 
step of processing a plurality of samples asyn- 

25 chronously through a plurality of measurement sta- 
tions so as to obtain at least one different physical 
property at each measurement station until each 
sample has been through each measurement sta- 
tion. 

30 7. A method of geophysical exploration com- 
prising the steps of: 

(a) obtaining a plurality of samples - of the 
earth's subterranean formations surrounding a well- 
bore at selected intervals along the length thereof; 

35 (b) processing the samples so as to obtain at 

least one physical property representative of their 
in-situ condition; and 

(c) producing an inverse log of at least one 
formation physical property representative of its in- 

40 situ condition along the length of the wellbore. 

8. A method of geophysical exploration, com- 
prising the steps of: 

(a) producing an inverse log of at least one 
selected physical property representative of its in- 

45 situ condition from samples of the earth's subterra- 
nean formations obtained at selected intervals 
along the length of the wellbore; 

(b) calibrating a set of geophysical survey 
data encompassing the wellbore location with the 

so inverse logs of the selected physical property; and 

(c) extrapolating the selected physical prop- 
erty to a location remote to the wellbore location 
with the calibrated geophysical survey data. 

9. The method of Claim 8 wherein samples are 
ss * obtained for each lithology traversed by the well- 
bore along its length. 

10. The method of Claim 8 wherein the sam- 
ples are obtained at intervals corresponding to the 
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at least the minimum resolution thickness of the 
geophysical survey data to be calibrated. 

. 11. The method of Claim 8 further including the 
step of smoothing the inverse log of a formation 
property with- a smoothing operator having a 5 
wavelength comparable to the minimum resolution 
thickness of the geophysical survey data to be 
calibrated. 
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